LANDAU MODEL FOR UNIAXIAL SYSTEMS WITH COMPLEX ORDER 
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We study the Landau model for uniaxial incommensurate-commensurate systems of the I class 
by keeping Umklapp terms of third and fourth order in the expansion of the free energy. It applies 
to systems in which the soft mode minimum lies between the corresponding commensurate wave 
numbers. The minimization of the Landau functional leads to the sine-Gordon equation with two 
r—{ • nonlinear terms, equivalent to the equation of motion for the well-known classical mechanical prob- 

lem of two mixing resonances. We calculate the average free energies for periodic, quasiperiodic and 
chaotic solutions of this equation, and show that in the regime of finite strengths of Umklapp terms 
only periodic solutions are absolute minima of the free energy, so that the phase diagram contains 
only commensurate configurations. The phase transitions between neighboring configurations are 
of the first order, and the wave number of ordering goes through harmless staircase with a finite 
^ ' number of steps. These results are the basis for the interpretation of phase diagrams for some ma- 

C/2 , terials from the I class of incommensurate-commensurate systems, in particular of those for A2BX4 

and BCCD compounds. Also, we argue that chaotic barriers which separate metastable periodic 
solutions represent an intrinsic mechanism for observed memory effects and thermal hystereses. 



64.70.Rh, 64.60.-i 



I. INTRODUCTION 



O ' Usual treatments of uniaxial incommensurate-commensurate (IC-C) phase transitions are based either on micro- 
scopic modelsjjsiith competing interactions or on phenomenological Landau theories. The relfiifant reviews can be 
found in Refs.EjD. The well-known example of the former is the Frenkel-Kontorova (FK) modeloo, iu which the wave 
T-H , number of ordering goes through the devil's staircase sequence of second order phase transitionsu. In the regime 
J> • of weak interactioBS the FK model can be continuated, and so reduced to the exactly solvable (i. e. integrable) 
sine-Gordon modeH. The solutions which then participate in the phase diagram are phase soliton lattices, i. e. com- 
mensurate regions separated by so called discommensurationso. The phase transition to the commensurate state is 
of the second (continuous) order, and the devil's staircase variation of the wave number is replaced by its simple 
■ continuous dependence on the control parameter. 
QQ The phenomenological Landau theory, another usual approach to the IC-C transitions, starts from the expansion 

^\ of the thermodynamic potential in terms of the order parameter, relied on the symmetry requirement by which the 
' order parameter is defined through one of the irreducible representations of the symmetry group of the normal ph. 
\ E. g., for structural phase transitions the order parameter is defined as a set of normal coordinates of the soft model 

Generally the minimum frequency of this soft mode may be located at an arbitrary point (i. e. star of wave vectors) 
I in the first Brillouin zone. The simplest irreducible representation for an uniaxial ordering is then two-dimensional. 
The corresponding basic ("minimal") form of the Landau expansion comprises, besides the leading normal terms, 
one, presumably the strongest, Umklapp term allowed by symmetry. This term favors a commensurate ordering and 
^ is responsible for the lock-in transition from the incommensurate ordering favored by the elastic term. Minimization 
;V ' of the Landau functioBjal again leads, after neglecting the space variations of the order parameter amplituddj, to the 
sine-Gordon equatiorfi-El, i. e. to the phase diagram equal to that of the FK model after the space continuation. 
The above approaches predict either a dense sequence of second order phase transitions (devil's staircase in the 
5— i _ FK model) or an isolated transition of the same tme. (Landau theory). Both possibilities are indeed close to the 



observations of IC-C transitions in some materialM^.^ A majority of materials however exhibits a more complex 
behavior comprising one or more first order phase transitions, memoriijeffects, wide ("global") hystereses, finite density 
of solitons at the very IC-C transition, etc (for a review see e. g. Ref.EJ). It is usually difficult to decide solely from the 
experimental observations, even for the most carefully prepared samples, whether such effects are of purely intrinsic 
or of some extrinsic origin. From the theoretical side, they cannot be explained within either of above approaches 
without extending the models. So far this problem was mainly considered by taking primarily into account some 
extrinsic agents, like external fields (e. g. electric field in ferroelectric materials), pinning centers, fixed or mobile 
defects, additional external periodic potentials with a periodicities different from that already present in the model), 
etc. 
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Another, more intricate possibility is that of intrinsic sources and mechanisms as the potential explanations for the 
afore mentioned phenomenally. In this respect the central question is the following; what are the simplest intrinsic 
extensions of the abovei-basic approaches which lead to phase diagrams with a finite sequence of first order transitions 
(i. e. harmless staircasellj) , and thus offer an inherent explanation for global hystereses and corresponding phenomena. 

The attempts in this direction were more successful in the realm of discrete moclels. The examples are models which 
include couplings-hetween next nearest neighbors, like .so-called DIFFOUR modelliJ, axial next-nearest neighbor Ising 
(ANNNI) mode]ll3'Ei.as well as vaxipus extensionstlTta, and models with two spin-like variables per site like those 
of Chen and WalkerEEl and JanssenEiJ. Both types of extensions were aimed mostly towards the interpretation of the 
phase diagrams observed in the family of A2BX4 compounds. 

On the other hand, the attempts within Landau models were based on the fotinai inclusiaruof more and more 
Umklapp terms (i. e. stars of wave vectors) into the basic models for the classes lE3~Eil and IlE3"E] of IC-C systems. 
From one side, the relevance of the Umklapp terms of high orders in the Landau expansion can be hardly justified on 
the physical grounds. Also, the ensuing analyses took into account only sinusoidal modulations, which, as the present 
study shows, is a too crude approximation for the determination of phase diagrams with harmless staircase, as well 
as for the interpretation of accompanying hysteretic effects. 

In contrast to such approaches, we propose in the present work a simple, physically well justified, extension of 
the basic Landau model for the class I, which is still framed within a "minimal" free energy expansion for a single 
star of wave vectors. The phase diagram which emerges from our model is characterized by a harmless staircase 
and first order transitions between highly nonsinusoidal configurations with different periods. Furthermore, a closer 
examination of configurations that participate in the phase diagram, and also of those that are not thermodynamically 
favored, enable a plausible explanation of the memory and hysteresis effects as the intrinsic (or at least semi-intrinsic) 
properties of IC-C systems. |— ■ 

Our considerations are based on a sine-Gordon model with two Umklapp termaSa. This type of model is physically 
well-grounded whenever the Landau expansion contains terms which favor two different commensurabilities, that are 
of comparable strengths. The most interesting case is realized with Umklapp terms of third and fourth order, tha 
lowest possible ones within the models with the Lifshitz invariant, appropriate for the so-called systems of the I classB 
(the systems of the II class have lock-in transition.s at the commensurabilities of order one and two, and are covered 
by essentially different type of Landau model^JcJ) . 

The mean-field (saddle point) approximation for our Landau functional leads to the Euler-Lagrange (EL) equation 
which has the form of the double sine-Gordon-Cfluation. This is one of the most intensively studied nonintegrable prob- 
lems in the contemporary classical mechanical tJ. The corresponding phase portrait contains periodic, quasiperiodic 
and chaotic trajectories, the latter appearing only when both nonlinear terms in the Landau functional are finite. 
As the strength of nonlinear terms increases, the chaotic trajectories occupy larger and larger portion of the phase 
space, destroying gradually quasiperiodic Kolmogorov-Arnold-Moser (KAM) layers, and eventually allowing only for 
some isolated periodic trajectories. The latter are orbitally unstable and therefore are not realized within the scope 
of classical mechanics. However, we show that just this tiny subset of the phase space comprises local minima of the 
free energy functional, i. e. the solutions (configurations) which participate in the thermodynappic, phase diagram. 
The question which then arises is analogous to that met in the analyses of the discrete modelfc'Eil, i. e. are there 
thermodynamically stable configurations among other, quasiperiodic and chaotic, trajectories. 

In order to analyze this additional, thermodynamic, aspect of the phase portrait, we calculate the average free 
energy for periodic, quasiperiodic and a representative set chaotic solutions of EL equation, with the aim to find, for 
given values of control parameters, those solutions that have the lowest value of the average free energy. We show that 
the chaotic configurations are never thermodynamically stable, in agreement with results obtained for some discrete 
modelso'Eil. The quasiperiodic configurations might be present in the phase diagram only when the Umklapp terms 
are weak enough, i. e. at temperatures slightly below the phase transition from the disordered to the incommensurate 
state. In the regime of strong Umklapp terms (to be specified later) the phase diagram is completely covered by 
periodic configurations, and the wave number of ordering passes through a finite number of values separated by the 
first order transitions, i. e. the corresponding staircase is harmless. 

The article is organized as follows. In Sec. II we introduce the Landau model of uniaxial ordering with two Umklapp 
terms and discuss its classical mechanical counterpart. The solutions of the Euler-Lagrange equation are considered in 
Sec. Ill, and the corresponding thermodynamic phase diagrams are presented in Sec. IV. Finally, in Sec. V we discuss 
possible implications to the phenomena observed in real materials, and compare our results with those obtained in 
the previous analyses of the similar models and other theories of uniaxial IC-C ordering. 
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II. MODEL 



We start from the assumption that the quadratic contribution to the Landau expansion has minima at wave numbers 
i+Q, —Q), where Q4 < Q < Qs, with Q4 = 27r/4 and Q3 = 27r/3. Here the unit length is taken equal to the lattice 
constant. The distances of Q from Q3 and Q4 are denoted by ^3 and 5^ respectively, with S3 + S4 = it/6 (Fig. |^). From 
now on we shall use S4 as an independent control parameter. Let us furthermore specify that the order parameter 
is complex, pe*"^. Limiting the further analysis to the temperature range well below the critical temperature for 
the transition from the disordered to the incommensurate phase, we also make the usual approximation of space 
independent amplitude pEl, and keep only the phase-dependent part of the free energy density. The latter reads 

fU^x) = h^f + Bcos\3(j) + 3(^-S4)x]+Ccos(4(j)-AS4x). (1) 
2 ax 6 



Here we scale the free energy functional 



F = y da:/(</),x), (2) 



by ^qP^, where ^0 is the correlation length in the x-direction. The first, gradient term in Eq. (|^) is the elastic 
contribution which favors the incommensurate sinusoidal ordering with the wave number Q. The second and third 
terms are the Umklapp contributions of the third and fourth order respectively. Due to the closeness of the wave 
number Q to both respective commensurate wave numbers, they are presumably the leading Umklapp contributions, 
provided both are allowed by symmetry. Their effective strengths are denoted by coefficients B and C which are 
proportional to the first and the second power of the amplitude p respectively. They are another two control parameters 
(beside 64) of the model (|l|). The temperature variation of p is expected to be the main source of the temperature 
dependence of B and C. 

The model (0) covers a variety of possibilities which may take place in particular physical examples. Besides the 
competition of each Umklapp term and the elastic term present already in basic sine-Gordon models, the essential 
new property of the present model is an additional competition between two Umklapp terms. The relative importance 
of these two terms is relied on both, the ratio of the strengths B and C and the position of the wave number Q, i. e. 
the ratio of S3 and S4, so that various regimes are possible. Regarding the expansion (|l|) it is reasonable to expect 
that the relative strength of two terms varies from the dominance of the third order term {B 3> C) at temperatures 
not far below the transition from the disordered phase to the comparable values of B and C at lower temperatures. 
However, even when B ^ C, the relative weakness of the fourth order Umklapp term can be compensated by a small 
value of 64 with respect to S3, i. e. by its much slower space dependence. In this case it is necessary to keep both 
Umklapp terms in the expansion (|l|). Although similar arguments may be invoked in favor of retaining some other 
pair of commensurate wave numbers, or even more than two of them, the example (|^) seems to be the most interesting 
one, due to the lowest possible powers of p present in the coefficients B and C. 

The configurations which take part in the thermodynamic phase diagram of the model (|l|) are the solutions of the 
Euler-Lagrange (EL) equation, 

0" + 3Bsm [30 + 3(- - S4)x] + 4Csin (40 - 4(54x) = 0, (3) 
6 

which for given values of the control parameters have the lowest value of the free energy averaged over the macroscopic 
length of the system, L, 

{F) = ^jdxmx),x]. (4) 

Before developing the appropriate method for the determination of such configurations, let us put few remarks on 
the equation (^. FxDiiiL.the classical mechanical side it represents the nonintegrable double resonance (i. e. double 
sine-Gordon) modeSo, with the corresponding Hamiltonian 

H{p^,(j),x) = ^-Bcos[30-f 3(^-(54)a;] - C cos (40 - 4(54a;) , (5) 

where = ^ = 0'. Obviously for i? = or C = Eqs. (^ and (|^) reduce to completely integrable sine-Gordon 
problems. For both B and C finite, one encounters the coexistence of two overlapping resonance domains. This can 
be easily seen with help of Poincare cross section. We introduce the auxiliary variable 
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(6) 



and plot the Poincare cross section in the phase space {ip,Pti,), ip = ipi^o + 3n),p^ = i/j'^xq + 3n). Here xq is the 
starting point of integration and n is an integer. The resonance domains are situated around eUiptic fixed points at 
(-0 = 0,p^ = ^m) and (0 = = j(2m+l)), where m is an integer. Their respective widths are 12VB/-k if C = 

and Yl\fC j-K if i? = 0. For small values of B and C (Fig. ^) the trajectories between two resonances conserve their 
topological form, while chaotic trajectories exist only very close to the separatrices of both resonances. As B and/or C 
increase (Fig. ^jb) the separatrices burst out into stochastic layers which grow and eventually merge between respnance 
domains. One gradually arrives at the threshold of the stochasticity (Fig. given by the Chirikov criterionEJ 

— (\/b + VC)«1, (7) 

at which the last KAM torus is destroyed, i. e. there are no more quasiperiodic solutions between two resonances. 
Chaotic trajectories are now free to diffuse through all phase space between two resonances. Let us here mention two 
points relevant for the further discussion. First, the widths of the chaotic layers grow exponentiallyBj as parameters 
B and C increase. Thus, the area between two resonances will be rapidly covered with chaatic layers. Second, KAM 
tori represent the main obstacles to diffusion of chaotic trajectories through phase space (E3 and references therein). 



III. SOLUTIONS OF THE EULER-LAGRANGE EQUATION 



Beside the classical mechanical context, our problem has an additional aspect, namely we are looking for the 
thermodynamically stable solutions, i. e. the trajectories in the phase space from Figs. |a,b:c which are local minima 
of the functional (^. Since we have to compare average free energies (^ of the trajectories present in the phase space, 
our first task is to specify numerical methods appropriate for the calculation of particular types of solutions. 

The orbitally unstable periodic solutions obviously cannot be determined by a direct integration of EL equations (|^), 
commonly used for calculation of orbitally stable solutions. It is therefore necessary to calculate them by using a 
suitable boundary value method for nonlinear equations. The most natural choice is the finite difference method, 
which is however rather demanding regarding the computer memory and time. It is therefore important to reduce the 
search for periodic solutions by establishing in an analytic way sufficient conditions on the possible values of periods. 
To this end we start from the relation 



</>(a; + F) = 0(a:)+(^P, (8) 

which holds for any periodic solution. Here P is its period and (/)p is the phase increment per period (note that the 
periodic solutions with finite (/>p belong to the rotational part of the phase space). Inserting Eq. (||) into the EL 
equation (||) taken at x + P one gets 

0" + 3P sin [30 + 3(- - 8Ax + 30p + 3(- - EAP\ + 4Csin (4(/i - ^b^x + 40p - ^b^P) = 0, (9) 
6 6 

where = 0(a;). Sufficient conditions on the values of parameters P and 0p follow from the requirement that Eqs. (|^) 
and (^) have the same form, i. e. that 

30p + 3(- -(54)P = 27rfc, 40P - 4(54P = -27rZ, (10) 
6 

where k and I are integers. Thus we get 

P = 4A: + 3;, ^p^b^P-l^. (11) 

Obviously, each periodic solution satisfciing the requirement ( [To| ) is uniquely defined by a pair of integers (A:, ^) which 
do not have a common integer factorcj. Note that the above procedure, in particular the step from Eq. (||) to 
the conditions ([lO|) , in principle does not forbid the existence of periodic solutions which do not belong to the set 
defined by Eqs. (|ll[). However, our attempts-|to locate numerically such solutions, although based on two independent 
algorithms, present and the alternative oneE3, always led to a negative result. This is an indication that the solutions 
with the periods ( |ll| ) are very probably the only possible periodic solutions, i. e. that the relations ( p^ ) are also the 
necessary conditions for their existence. 

The solution (^(x) with the period (O) has the total wave number (that measured from the origin of Brillouin zone) 
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(12) 



The values of q allowed by conditions ( [lO| ) form a Farey tree, shown in Fig. |^ for the wave numbers between q — 1/3 
{k = 0, I — 1 and P — 3) and q = 1/4 {k = 1, I = and P = 4). Thus, already at this introductory stage of the 
analysis we conclude that the model (^ has the phase -d-iagram with branchings between neighboring commensurate 
configurations equivalent to those of ANNNI modelsli3'Ea. 

The periodic solutions q = 1/3 {k ^ 0, I ~ 1 and P — 3) and g = 1/4 (fc = 1, ^ = and P = 4) in the Farey 
tree of Fig. ^ are the basic commensurate configurations, belonging to the Umklapp terms of third and fourth order 
respectively. The wave numbers at lower levels of the Farey tree from Fig. |^ represent higher order commensurate 
solutions which correspond to all positive values of k and I. They are situated between two main resonances in the 
Poincare cross section shown in Fig. ^. Note that for small values of B and C their positions in the phase space 
(Fig. ||a) perfectly match positions in the Farey tree (Fig. H). As the parameters B and C further increase, the 
positions of periodic solutions that are embedded in chaotic layers become slightly intermixed since there are no 
more KAM tori between two resonances that restrict their positions in phase space (Fig. |2|3). We do not include the 
parts of the Farey tree belonging to negative values of k and/or I since, as it will become clear later, they are not 
thermodynamically stable for < 84 < tt/G and B,C > 0. 

In the next step we specify boundary conditions for a particular periodic solution (j)ki{x). Since every periodic 
solution possesses at least two inflection points, we chose one of them, a;o, to be the initial point of integration, i. e. 
the left end point of one period. Thus 4>'ki{x — x^) = 0. The boundary conditions now read 

(l)ki{x ^ xo + P) = (f>ki{x = Xo) + 0P, 

0h(x Xo + P) = (j!)fc,(x = Xo). (13) 



Since the values of of P and (pp follow from the choice of integers {k,l) (Eqs. (|11|)), it remains to establish the 
connection between the other three parameters, xo, (pkiix = xq) and (f>'f^i{x = xq), which figure in the conditions (p^). 
As it follows from the EL equation (^) with x = xq, 

3B sin [30(xo) + 3(^ - S4)xo] + 4Csin (4(?i(xo) - 464X0) = 0, (14) 

Xo and (/'(xo) are not independent. Even more, the numerical experience suggests that for a given periodic solution 



both Xo and (/)(xo) do not vary as we change B and/or C, i. e. that Eq. (14) in fact decomposes into two conditions 



TT 

3(/>(xo) + 3(- - S4)xo = Mix, 
6 

4(/)(xo) - 464X0 = -Ntt, (15) 
where M and N are integers. This means that xo and '/>(xo) may have values 

xo = i(4M + 37V) (16) 

and 

0(xo) = 54x0 - JiV. (17) 

These relations would allow for 2P values of Xo and an infinite number of values for (f>{xo) (for a general value 
of ^4). The further analysis of symmetry properties of the problem (^), as well as the numerical insight, however 
indicate that for any choice of periods P and (pp this enumerable set is highly degenerate and reduces to only two 
distinct (nondegenerate) solutions. The convenient choices of xq and (j){xo) characterizing these solutions for various 
combinations of odd and/or even values of integers k and I are listed in Tab. ^ 

The above analysis simplifies drastically the numerical procedure, since after specifying the parameters k, I, xq and 
<j){xo), the determination of a given periodic solution follows from the variation of the single remaining parameter 
(j)'{xo)- In accomplishing this procedure it appears convenient to eliminate, by transforming the variable <j){x), the 
explicit X— dependence from one of Umklapp terms in the EL equation (|^), and to keep this dependence in the term 
with a weaker amplitude. Thus for B larger than C we use the variable ip{x) (Eq. (^)) instead of (j){x), so that the 
EL equation 

iP" + 3B sin (3?A) + 4C sin (4^/^ - — x) = 0, (18) 
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and its solutions ip{x) do not depend on ^4. The corresponding free energy then acquires a 54-dependent term in the 
form of Lifshitz invariant, 



I dx{^[ij' -i^-Si)]^ + Bcos{3^) + Ccos{Ai;-^x)}, (19) 

which simphfies the calculation of the (54-dependence of the averaged free energy for any particular periodic solution 
of the EL equation. As is visible in Fig. the form of periodic solutions resembles to that of multisoliton solutions 
of simple sine-Gordon model (still, note the slight modulation of commensurate regions, i. e. between discommen- 
surations) When C is larger than B, it is appropriate to introduce an analogous variable which makes the C-term 
X— independent, namely x{x) — 4>{x) — S^x. Again, the corresponding EL equation does not depend on ^4. The 
boundary conditions have to be modified correspondingly for both transformations. 

Although the steps described above greatly simplify the numerical procedure, the finite difference method poses the 
limitations on the computer memory and time which do not allow us to calculate solutions with periods well above 
100. Note in this respect that the nonlinearity of EL equation (|^) or ( [T^ ) forces us to use about 1000 mesh points per 
period in order to get solutions which are reliable enough. 

Periodic solutions of EL equation ^ show several interesting properties which are important for analysis of phase 
diagrams. We notice that for some values (or ranges of values) of parameters B and C one of the two periodic solutions 
with the same values of k and I from Tab. || ceases to exist (see Fig. |]). In general the solutions with the lower value 
of averaged free energy are more robust to this disappearance. We do not go into a closer analysis of this effect, but 
only indicate that it seems to be closely connected with the destruction of KAM tori as B and C increase. 

Another interesting property of periodic solutions is the splitting in averaged free energies of two solutions with the 
same values of (fc, /) (see Fig. H). As the parameter C gradually increases from zero, while keeping B fixed, values of the 
difference between these two energies increases, thus making one periodic solution more and more thermodynamically 
favorable with respect to the other. This splitting is larger for the solutions with smaller periods. The qualitative 
consequence is that such solutions participate over greater and greater parts of the phase diagram as parameters B 
and C increase. 

For the calculation of quasiperiodic and chaotic trajectories we use standard, Adams or Runge-Kutta-Merson, 
metho^ for an initial value problem. Quasiperiodic trajectories, as a building objects of KAM tori, are orbitally 
stabler. Chaotic orbits, although certainly orbitally unstable, are diffusive through all the corresponding chaotic 
layer in the phase space, so that by picking one of them we get practically the averaged free energy for all chaotic 
solutions in that layer. Thus, in order to calculate the averaged free energies of quasiperiodic and chaotic solutions 
we chose initial values by random (probability of picking periodic solution instead of quasiperiodic or chaotic ones is 
equal to zero), and carry out the integration as long as the accuracy is satisfying. The fact that the averaged free 
energy of quasiperiodic and chaotic solutions can be determined only to a limited accuracy was already pointed out 
by Fradkin et. ala who estimated the degree of accuracy for a given type of solution. i_. 

The estimation of the common averaged free energy of chaotic solutions within a given layer can be done as foUowsEj. 
The average value of Umklapp terms in the expression ( pj| ) is zero since these terms contain trigonometric functions 
with an argument which chaotically (randomly) varies with x. For the fourth order Umklapp term cos {4:^ — ^x) we 
have 

27r 2n 
(cos {—x) cos (4V')) = (sin {—x) sin (4^^)) = 0, (20) 
o o 

while for the third order Umklapp term we have an average of cos which is also zero. The averaged free energy is 
thus given by the integral of the gradient term ^[ip' — (tt/G — ^4)]^. The latter depends on the position and the width 
of chaotic layer in the phase space, i. e. only on the dependence of ip' on a; along the trajectory in this layer. 

In order to determine a solution with the lowest average free energy we follow solutions (periodic, quasiperiodic 
and chaotic) with initial conditions that belong to the line {iIj{xo) = 0, ip'{xo)) in the phase space {ip, ip') (Fig. ||), and 
compute their average free energies. For small values of B and C, periodic and quasiperiodic solutions are regularly 
arranged in the phase space (Fig. ||a), with hardly distinguishable average free energies (Fig. ^). In order to show 
that the solution with the lowest free energy is periodic, we follow downwards the branch of the Farey tree (Fig. |^) 
which starts at the point with the averaged free energy lower that those for neighboring points above and below this 
point. It is a numerical evidence that the average free energies increase (and tend to some finite value) as we go down 
through successive branch points, i. e. through the solution with larger and larger periods. Quasiperiodic solutions 
can be regarded as asymptotic limits of series of periodic solutions defined by successive branchings in the Farey tree, 
in which the period and the phase increment tend to infinity (but with a finite, irrational, value of q). The averaged 
free energies of quasiperiodic solutions thus should be equal to the limiting values of averaged free energies at a given 
branch, which are, as is argued above, higher than the averaged free energy of the starting periodic solution. Since 
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this argument is based on numerical calculations, it cannot be extended to very small values of B and C for which 
the solution with the lowest free energy, as well as the solutions at the accompanying branch in the Farey tree, have 
too large periods. 

In the range of intermediate values of B and C (Fig. ^3) there are intervals of initial conditions in which quasiperiodic 
solutions disappear, and only chaotic and periodic solutions are present. The chaotic layers can be easily recognized 
in the Poincare cross section (Fig. |^). The average free energies of periodic solutions then look as needle-like minima 
immersed in the average free energy of chaotic layers, represented by plateaux in Fig. Finally, for large values 
of B and C (Fig. |6|3) for which the Chirikov criterion (Q) is fulfilled, there remains a single chaotic layer between 
two resonance domains (Fig. |^c), while the number of existing periodic solutions gradually decreases as B and C 
increase. Since there remains a finite number of corresponding well defined needle-like minima, it is sufficient to limit 
the numerical calculations to the search for existing periodic solutions, and to find out among them the solution that 
has the lowest average free energy. 



IV. PHASE DIAGRAM 



We have argued in the previous section that the configurations with minimal average free energy are among periodic 
solutions of EL equation (^ . Before presenting results of numerical calculations which confirm this expectation, we 
briefly discuss the parameters present in the model (^. 

The parameters B and C depend on external conditions, most usually on temperature and pressure. As it was 
mentioned in Sec. II, they depend on the amplitude of the order parameter linearly and quadratically respectively. At 
temperatures closely below T/, the temperature of phase transition from the disordered to the incommensurate phase, 
the ratio B/C is proportional to (T/ — T)^^/^. A more complete insight into the temperature dependence of the order 
parameter, and of the ratio B/C as well, in the wider temperature range below T/, can berebtained from the neutron 
scattering, NMR and similar experimental data for particular materials (e. g. referenceaSafij) . As for the pressure 
dependence of B and C, it can be specified only after the insight into the microscopic model for a particular material 
on which the Landau theory is based. The parameter 64 also might be temperature and/or pressure dependent. 
Usually, in a concrete physical situation certain dependences may be regarded as dominant. E. g., when temperature 
varies and pressure is constant ^4 can be often regarded as constant, while B and C are temperature dependent. 
Having this in mind we simplify the further discussion by keeping one of the parameters fixed and concentrating on 
phase diagrams in remaining two-dimensional parameter subspaces. 

The role of the parameter 64, the position of the instability with respect to the wave number of the fourth order 
commensurability, is expressed through the Lifshitz invariant d4'4>'{x) in Eq. (19) which favors the incommensurate 
ordering. On the other side two Umklapp terms favor commensurate orderings with their respective wave numbers. For 
64 —^ 0, and fixed values of parameters B and C, the Umklapp term of the fourth order dominates with respect to that 
of the third order, and the thermodynamically stable periodic solution is expected to have the wave number qq — 1/4. 
On the same footing, for ^4 near it/6 (i. e. for ^3 0) the stabilization of the modulation with qq — l/i is preferred. 
For < 64 < tt/6 we expect that some other higher-order wave numbers of modulation become thermodynamically 
stable and that they follow the order specified by the Farey tree from Fig. ||. 

Let us now fix the parameter B and allow for the variation of the parameters 64 and C. For a particular value of C 
we find periodic solutions of the EL equation (|l^) by following the steps from Sec. Ill, and calculate their average free 
energy ( p^ ) for a relevant range of values of the parameter 64. Then we determine a solution which is the absolute 
minimum of the average free energy for a given value of ^4, and in particular the isolated values of 64 at which first 
order phase transitions take place since two (or more) configurations are simultaneously absolute minima of the free 
energy. Varying also systematically the parameter C we obtain the phase diagram, as shown in Fig. ^ for B = 0.02. 
All lines in this diagram represent the phase transitions of the first order between the periodic configurations with 
different wave numbers (which are denoted only for few dominant phases in the diagram). Note that the Chirikov 
line (|^ is at C « 0.0145, and that below C « 0.01 there is a proliferation of configurations with commensurabilies 
of higher and higher orders. The absence of these configurations at larger values of C is mostly due to the fact that, 
although they exist as solutions of the EL equation, their average free energies are too high in comparison to those 
of the solutions with lower commensurabilities. In addition, some periodic solutions simply cease to exist as C (or 
B) increase, as shown in Fig. ||. Note also that only one of two different classes of periodic solutions with the same 
values of k and I participates in the phase diagram in Fig. i. e. that characterized by the initial conditions from 
the second rows (depending on evenness and oddness of integers k and I) in Tab. ^. Still, we find out numerically that 
the average free energies for two different solutions with the same (fc, I) may change order, i. e. that the solutions 
from the first rows in Tab. | may have lower average free energy than those from the second rows provided they are 
of rather high commensurability. Thus, it is somewhat surprising that in the phase diagram in Fig. ^ the periodic 
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solutions with only one type of boundary conditions prevail. We shall come back to this point later in Sec. V. 

In addition to the phase diagram, we plot in Fig. ^ the corresponding staircase, i. e. the wave number of the stable 
configuration vs parameters C and S^. As long as C is not very sniall there is a finite number of steps, i. e. we obtain 
the so-called harmless staircase, introduced by Villain and Gordonll3. We stress that the most interesting property of 
the phase diagram from Figs. |^ and ||, the presence of a finite (small) number of stable commensurate configurations, 
is encountered in the regime of rather high values of parameters B and C. The phase portrait of the EL equation (j^) 
is then almost everywhere chaotic (Fig. 0c) and there are no more quasiperiodic solutions between two resonances. 
By increasing further the values of parameters B and C one eventually comes to the phase diagram in which only two 
main commensurate phases {q = 1/3 and q = 1/4) take place. 

Another possible presentation of the phase diagram is that with a fixed value of the parameter and with varying 
parameters B and C. It is presumably closer to usual physical situations in which only weak temperature and 
pressure dependences of 64 are expected. The construction of the (_B, C) phase diagram is however computationally 
more demanding, since one has to look for the solution with the lowest average free energy within a set of solutions 
for given values of B and C, i. e. one has to calculate the whole set of periodic solutions of the EL equation [Eq. ^ 
or (y_8|)] for each point in the two-dimensional phase diagram. To this end we use a mesh of points which is dense 
enough in the {B, C) plane, and determine the solution with the lowest average free energy at each point. The phase 
diagram obtained in this way for ^4 = 7r/12 is shown in Fig. ^. Note that again only configurations with rather low 
orders of commensurability, i. e. with small values of parameters (k, I), are present above the Chirikov line (Eq. (0)), 
while below this line the diagram is more complex since a great number of first order transitions takes place within a 
small part of the phase diagram. 



V. CONCLUSION 



The most important conclusions of the above analysis follow from the thermodynamic phase diagram obtained in the 
regime of comparable strengths of two Umklapp terms included into the Landau expansion (|^) . At first we emphasize 
that only one type of solutions of the corresponding EL equation, namely periodic configurations, participates in the 
phase diagram. Furthermore, all phase transitions between successive commensurate phases are of the first order, 
so that the wave number of ordering follows a harmless staircase with a finite number of steps. The examples of 
such phase diagram, namely series of successive iock-in commensurate-commensurate transitions with accompanying 
effects which characterize first order transitionstj, are often encountered in particular materials. Here we focus our 
attention on few well known examples. 

One of the most studied type of materials are A2BX*-<;oKipounds, among which we take Rb2ZnBr4 as a prominent 
representative. Early neutron diffraction measurementsE^ES of the temperature variation of modulation wave number 
revealed the existence of several higher-order commensurate phases. The more complete pressure-temperature phase 
diagram followed., from various subsequent data, in particular again from the neutron diffraction measurements of 
Parlinski et al.EJa. It resembles to a great extent to our phase diagrams from Figs. ^ and ^. Note also that 
the phenomenological formula for wave numbers of observed commensurate phases introduced by Parlinski et al.c^l 
coincides to our expression for Farey tree (|l2|), which is, as is shown in Sec. Ill, inherent to the-|model (|^). Harmless 
staircase are clearly seen in e. g. pressure variation of the wave vector for a fixed temperaturec3, with steps going as 
1/3, 7/24, 2/7 and 1/4 by increasing pressure. They are accompanied by hysteresis in pressure and temperature runs, 
which are particularly strong when only a few steps appear in the phase diagram. This corresponds to the regime of 
rather high values of parameters B and C, in which the phase diagram contains only few commensurate phases and 
the average free energy of chaotic plateau is well above average free energies of periodic solutions (Fig. ^). 

Existing theoretical approaches to the (in) commensurate orderings in A2BX4 compounds, in particular ta_jthe 
appearance of series of commensurate phases, are most |ly . . phen omenological , based either on Landau expansionsCJ or 
on the discrete models of competing local interactions" ^r^rn. The justification for the continuous Landau models, 
which are generally appropriate to-weak coupling systems, comes from many experimental indications, starting from 
the early neutron scattering dataKlEa, showing a well-defined dispersion curves for collective modes with distinct 
soft-mode minima. However, the previous analyses of Landau models were restricted to purely sinusoidal modulation, 
and, as such were not able to explain the appearance of phases with commensurabilities of orders higher than three 
or four. It was therefore prop|Osed that such phases appear due to the presence of Umklapp terms of higher orders 
in the free energy expansionE3~Ej. This explanation, which is based on the assumption that distinct commensurate 
stars of-Wiaye vectors are necessary for the stabilization of, presumably sinusoidal, phases with corresponding wave 
vectorscS'LJ, is not convincing since the Umklapp terms of order higher than four are expected to be negligible in 
weakly coupled systems with a displacive order. 

For these reasons the more recent attempts turned again towards another type of approaches, those which assume 
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strong couplings, so that the lattice discreteness has to be taken into account. Originally the sequences of IC-C and 
C-C phase transitions were within this scheme interpreted in terms of the FK model as devil's staircase dependences of 
the wave number of ordering, i. e. as dense sequences of second order phase transitionsH. However, observed staircase 
rarely resemble, even within experimental limitations, to the dense devil's staircase. Beside, phase transitions between 
successive C phases are usually of the first order. The phase diagrams which are closer to experimental findings 
may be bawever obtained by various extensions of the basic FK model, e. g. by includiag an additiopalJjarmonic 
potentialuScZl. Also, the more complex models of competing interactjons, e. g. DIFFOURfi^ and ANNNIHTE^ models, 
as well as models that assume two critical modes per lattice siteEj^Ej, are particulady successful in describing the 
phase transitions in A2BX4 compounds. Within some of these models (e. g. Refs.c2l andcj) one also obtains the 
first order phase transitions between configurations having the same wave numbers but different symmetries. As was 
already stated in Sec. IV, this is not the case within the model (l|), i. e. although the EL equation (||) may possess 
two types of solutions with the same periods, only one type of solutions participates in the phase diagram. 

The present analysis again starts from the minimal Landau expansion (with terms up to the fourth order), but 
takes into-consideration all solutions of the corresponding EL equation. In particular it indicates that the theoretical 
approaches, proposed together with the first neutron scattering measurements on Rb2ZnBr4, might be essentially 
sufficient for the understanding of complex phase diagrams in A2BX4 materials. The more detailed analysis which 
takes into account some additional aspects, like the couplings to the homogeneous polarization and strain which 
appear in some materials as secondary order parameters, will be done elsewhere. 

Betaine-calciumchloride-dihydrate (BCCD), together with its deuterated version D-BCCD, belongs to the second 
type of intensely studied materials with the commensurate lock-ins. It shows_afi_exceptionally rich staircase going 
from (7 = 1/3 down to 9 = 0, with numerous intermediate steps of higher orderal3~Ej. A closer insight into the region 
of phase diagram with the wave number between g = 1/3 and g = 1/4 shows that only upper right triangle of the 
Farey tree from Fig. |^ is realized, i. e. the phase diagram is mostly covered by wave vectors close to q = 1/4, and 
not by those close to q = 1/3. This sequence of IG-|P-|transitions was successfully interpreted within various discrete 
models with competing interactions, e. g. in Refs.EllE3. Within the model (|l|) such phase diagram corresponds to the 
regime in which the fourth order Umklapp term dominates with respect to that of the third order. Also, two types of 
extensions of our model may lead to the stabilization of commensurate phases with q < 1/4. Namely, one may allow 
for negative values of the parameter S4, or start with other Umklapp terms, e. g. with those of the fourth and fifth 
order, and pursue the analysis analogous to that of Sees. Ill and IV. 

We also mention some other materials that exhibit a sequence of IC-C and C-C phase transitions between q = 1/3 
and q = 1/4, but are not so extensively studied as previous two examples. E. g., Denoyer et alE3 investigated 
NIl4lISe04 and its deuterated version ND4DSe04 by neutron diffraction, and found the harmless staircase and first 
order phase transitions, accompanied by the coexistence of several phases in the relatively wide range of temperatures. 
A series of IC-C and C-C phase transitions are also observed in BaZnGe04 its. X-ray diffraction measurements by 
Sakashita et al.Eil, and in electron diffraction measurements by Yamamoto et al.E3 which also provide dark field images 
of discommensurations appearing in the vicinity of q = 1/3 phase. An example of particularly sharp transition from 
q = 1/3 to g = |lY4, with a very wide temperature range of the coexistence of these two commensurate phases, was 



ongs to A2BX4 family. 

who also studied continuum systems with competing 



found by BrodaE3 in (NIl4)2CoCl4, the material that also 

The free energy (|l|) is similar to that of Fradkin et al.l 
periodicities. The only difference between two expressions is the absence of the factors 3 and 4 in front-of the variable 
tf>{x) in the cosine terms of the model^j. However, in contrast to ours, the analjsis carried out in Ref.O is limited to 
the close vicinity of the separatrices (and hyperbolic points) in the phase spaceEJ, i. e. to the dilute soliton lattices. 
Then thepContinuum model can be converted into a discrete mapping of the FK type, analyzed in detail previously 
by Aubryd. Our analysis covers the whole phase space, i. e. all solutions of the EL equation (^, and in particular 
the whole class of periodic-configurations. In particular, our thermodynamic phase diagram (Figs. ^ and |^) includes, 
in contrast to that of Ref.u, the most interesting part of the phase space, namely that between two resonances (i. e. 
sets of hyperbohc points). _ 

The modelEj was the starting-^oint in the investigationEll of the memory effects in systems with IC modulations, 
based on the earlier propositionCi that mobile defects might be responsible, by forming defect density waves, for the 
sensitivity of the IC ordering on the thermal history of crystal, observed e. g. in thioureaEHl. ErrandoneaEll argued 
that the double sine-Gordon model, with two lock-in potentials originating from the lattice the defect density wave, 
is an appropriate description of this phenomena. 

The model (|l|) propiies the explanation of memory effects (together with thermal hystereses), without referring, 
in contrast to modelstZrEj, to defects as an intrinsic ingredient of the theory. At first, we note that the crossings pf 
lines of first order phase transitions in Figs. 0, ^ and |^ are accompanied by hystereses. Our preliminary analysisE3 
indicates that these hystereses may be rather wide on e. g. temperature scale. Furthermore, the present analysis 
of the EL equation (||) shows that periodic solutions, which constitute the phase diagram, are immersed as isolated 
points into the environment of chaotic configurations. This environment prevents both, the continuous variation of 
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the wave number of ordering and the continuous phase transitions of the second and higher orders. The average free 
energy of chaotic solutions from Fig. ^jc is the measure of the energetic barrier which the system has to overcome in 
order to pass from some periodic (metastable) configuration to another one with lower free energy. This is expected 
to be a common property of models which are nonintegrable (beside being nonlinear), and have thermodynamically 
stable periodic configurations isolated in the chaotic phase spaceEp. 

The, memory effects are also observed in class II of IC systemstB. The detailed analysis of phase diagram for this 
clasg23 led to the conclusion that the corresponding phenomena seen in particular materials may be as well interpreted 
in terms similar to those presented above. However, it was also stressed lB that defects may have a secondary role 
as triggers which favor the stabilization of some domain patterns. This interpretation invokes neither the mobility of 
defects nor the formation of defect density waves. The analogous secondary influence of defects on memory phenomena 
is expected also in presently investigated systems of class I. 

Finally, let us mention a conuaon problem that arises in the analysis of continuous nonintegrable Landau models 
for uniaxial systems of classes ]n3 and IIc2l, in which periodic solutions have an essential role in the extremalization 
of corresponding thermodynamic functionals. We remind that there is no firm universal principle which would favor 
the thermodynamic stability of the (meta)stable periodic configurations on account of other, quasiperiodic or chaotic, 
solutions of EL equations. Some hints in this direction for "autonomous" functionals (those-.for which the free energy 
density does not depend explicitly on x) follow from the recently derived general criteriaEj based on the additional 
extremalizations (like e. g. those involving boundary conditions). However, these criteria cannot be directly applied 
to the present nicijdel since the explicit x-dependence in Eq. (^ introduces fundamental singularities in the additional 
extremalizatonsEHl. Thus, the most important property of the phase diagrams from Figs, [t] and ||, their complete 
coverage by a finite number of periodic configurations, still awaits for a deeper understanding. 

We acknowledge numerous useful discussions with V. Dananic. The work is supported by the Ministry of Science 
and Technology of the Republic of Croatia through the project no. 119201. 
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FIG. 1. Brillouin zone with the soft-mode minimum at Q, and the commensurate wave numbers of third (Q3) and fourth 
(Q4) order. 



FIG. 2. The Poincare cross sections for the Euler-Lagrange equation ( |l8| ) and the choice of parameters: xq = 1.5, ip{xo) — 0, 
B = C = 0.002 (a), B = 0.008, C = 0.006 (b), B = C = 0.02 (c). The period which defines the section is equal to 3. Symbols 
for the periodic solutions {k,l) are: sohd Q (0,1), solid □ (1,0), solid O (1,1), solid A (1,2), solid < (1,3), solid V (1)4), solid 
t> (2,1), + (2,3), X (3,1), * (3,2). 



FIG. 3. Farey tree for wave numbers q defined by Eq. (|l3). 



FIG. 4. The periodic solutions il){x) from the class (1,/). The parameters are: B = 0.02, C = 0.02, xo ~ 1.5, tp{xo) = 0. 

FIG. 5. Average free energies of the periodic solutions from the class (1, 1) as the function of C for B — 0.02 and ^4 = 7r/12. 
Figure (b) is the enlarged detail of the figure (a) with energies lower than 0.01. Solutions with lower average free energies are 
those from the second rows in Tab. ^ Note from the figure (b) that e. g. the upper solution (1,6) does not exist for few 
subranges of the values of parameter C, and that both solutions from this class cease to exist for C > 0.03. 

FIG. 6. Average free energy vs tp'o of periodic (A), quasiperiodic (Q) and chaotic (solid O) solutions, for B — 0.002, C — 0.002 
(a), B = 0.008, C = 0.006 (b), B = 0.02, C = 0.02 (c), and 84, = 7r/12, xq = 1.5. The {k,l) indices for the periodic solution 
with the lowest average free energy (solid A) are (3,4) in the figures (a) and (b), and (1, 1) in the figure (c). The insets in 
figures (a) and (b) are enlarged neighborhoods of the free energy minima. 



FIG. 7. Phase diagram in the (C, 84) plane for B — 0.02. The numbers in the figure are periods of some stable commensurate 
phases. The dashed line at C ~ 0.0145 represents the Chirikov criterion (0). 

FIG. 8. The wave number of modulation qo vs C and ^4 for B — 0.02. The dotted cross-section represents the Chirikov 
criterion (^. 

FIG. 9. Phase diagram in the {B, C) plane for ^4 = 7r/12. The numbers in the figure are periods of some stable commensurate 
phases. The dashed curve represents the Chirikov criterion n^. 
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